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Abstract 



Using a group theoretical approach we derive an equation of motion for a mixed 
quantum-classical system. The quantum-classical bracket entering the equation pre- 
serves the Lie algebra structure of quantum and classical mechanics: The bracket is 
antisymmetric and satisfies the Jacobi identity, and, therefore, leads to a natural de- 
scription of interaction between quantum and classical degrees of freedom. We apply 
the formalism to coupled quantum and classical oscillators and show how various ap- 
proximations, such as the mean-field and the multiconfiguration mean-field approaches, 
can be obtained from the quantum-classical equation of motion. 
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1 Introduction 

Many phenomena in nature are described by quantum mechanics at a fundamental level and 
with high precision. Yet, there exist numerous situations where mixed quantum-classical 
models are needed. In some cases the phenomena are too complex to allow for a fully quan- 
tum approach, in others a consistent quantum theory is lacking. Classical mechanics often 
provides a more suggestive description and a clearer picture of physical events. Applica- 
tions of various quantum-classical approaches range from biochemical and condensed matter 
chemical reactions, where the large dimensionality of the systems of interest requires ap- 
proximations, to the evolution of the universe and cosmology, where no theory of quantum 
gravity has been established. 

The issue of treating quantum and classical degrees of freedom within the same formalism 
was recently discussed in a number of publications 0, Q . The interest was spurred by the 
cosmological problem of defining the backreaction of quantum matter fields on the classical 
space-time background, where classical variables should be independently correlated with 
each individual quantum state. The traditional approach fails to satisfy the last requirement. 
(For a fully quantum approach to cosmology see reference [Q.) Somewhat earlier a similar 
situation was encountered in chemical physics, where quantum-classical trajectory methods 
were employed to model gas phase scattering phenomena and chemical dynamics in liquids |^ 
^, 1^, H, ^, [10[ . It was noticed in these studies that asymptotically distinct quantum evolutions 
should correlate with different classical trajectories. 

The first relationship between quantum and classical variables is due to Ehrenfest who 
showed that the equation of motion for the average values of quantum observables coincides 
with the corresponding classical expression. (Surprisingly, the first mathematically rigorous 



treatment on the subject was not carried out until 1974, see reference [12|.) Ehrenfest's result 
leads to the mean-field approach, where classical dynamics is coupled to the evolution of the 
expectation values of quantum variables [|T3], |T^, [T^, The mean-field equations of motion 
possess all of the properties of the purely classical equations and are exact as far as the 
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mean values of quantum operators are concerned. However, an expectation value does not 
provide information of the outcome of an individual process. The mean-field approach gives a 
satisfactory description of the classical subsystem as long as changes within the quantum part 
are fast compared to the characteristic classical time-scale. If classical trajectories depend 
strongly on a particular realization of the quantum evolution, the mean field approximation 
is inadequate. The problem can be corrected, for instance by introduction of stochastic 
quantum hops between preferred basis states with probabilities determined by the usual 
quantum-mechanical rules |]T^ 0, P^ . 

Similarity between the algebraic structures underlying quantum and classical mechanics 
provides a consistent way of improving upon the mean-field approximation, as explored in 
references 0, |5, • In those studies the aim was to derive a quantum-classical bracket that 
reduces to the quantum commutator and the Poisson bracket in the purely quantum and 
classical cases. In addition to the reduction property the bracket should satisfy other criteria 
so as to give physically meaningful pictures of quantum-classical evolutions. In particular, an 
antisymmetric bracket conserves the total energy and a bracket satisfying the Jacobi identity 
ensures that the Heisenberg uncertainty principle is not violated. 

Recently, one of us (VVK) proposed a natural mathematical construction, which 
we name p-mechanics, enveloping classical and quantum mechanics. Formulated within 
the framework of operator algebras, the p-mechanical equation of motion reduces to the 
appropriate quantum or classical equations under suitable representations of the algebra 
of observables. In this paper we extend the ideas of p-mechanics to incorporate mixed 
quantum-classical descriptions. In particular, we derive the quantum-classical bracket and 
explicitly show that it satisfies the properties common to quantum and classical mechanics. 
Using the technique described it is possible to construct families of mixed quantum-classical 
approaches, each having a specific set of properties. 

The format of this paper is as follows: In Section ^ we summarize p-mechanics and 
introduce the essential mathematical definitions. In Section ^ we construct the simplest 
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p-mechanical model that adopts two distinct sets of variables associated with quantum 
and classical degrees of freedom. By taking the appropriate representation we derive the 
quantum-classical bracket and show that it is antisymmetric and obeys the Jacobi identity, 
that is, it possesses the two major properties shared by the quantum and classical brackets. 
In Section ^ we work out the case of coupled classical and quantum harmonic oscillators. 
Finally, in Section ^ we discuss how various approximations to the general quantum-classical 
description can be obtained, including the mean-field and the multiconfiguration mean-field 
approaches. 

2 P-mechanics 

2.1 The Elements of P-mechanics 



We recall the constructions from references EQ, 21] together with appropriate modifications 



Definition 1 An operator algebra ^ gives a p-mechanical description of a system if 
the following conditions hold. 

1.1. The set ^ of all irreducible representations iTh of ^ is a disjoint union of subsets 
^ = Upgp^p parameterized by the elements of a set P. The elements of the set P 
are associated with different values for the Planck constant. We refer to this set as 
the set of Planck constants. If for po the set consists of only commutative (and, 
therefore, one-dimensional) representations, then ^Pp^ gives a classical description. If 
^po — {^po} consists of a single non- commutative representation tt^q, then ^Pp,, gives 
a purely quantum model. Sets ^p of other types provide mixed (quantum-classical) 
descriptions. 

1.2. Let ^ be equipped with a natural operator topology (for example, it may be the 
Jacobson topology or the *-bundle topology ill). Then P has a natural factor 
topology induced by the partition *P = Upgp. 
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1.3. (Dynamics) The algebra ^ is equipped with the one-parameter semigroup of transfor- 
mations G{t) : ^ ^ t G M"*". All sets p G P are preserved by G{t). Namely, for 
any vr G all new representations ttj = vr o G{t) again belong to ^p. 

1.4. (The Correspondence Principle) Let S : p ^ S{p) G be an operator- valued section 
continuous in the *-bundle topology [E^, |23 over P. Then for any t, i.e., at any 



moment of time the image St{p) = G{t)S{p) is also a section due to statement |L3|. In 
the *-bundle topology the sections St{p) are continuous for all t. 

Having listed these quite natural conditions, we do not yet know how to construct p- 
mechanics. Next, we describe an important particular case of group quantization [^]. All 
components of p-mechanics (operator algebra, partition of representations, topology) readily 
arise there. 

Construction 2 Group quantization consists of the following steps. 

2.1. Let Q = {xj}, 1 < j < N he a set of physical variables defining the state of a classical 
system. Classical observables are real-valued functions on the states. 

The best known and the most important case is the set {xj = qj,Xj^n = Pj}^ 1 ^ J ^ 
n, N = 2n of coordinates and momenta of classical particles forming an n degree of 
freedom system. The observables are real valued functions on M^*^. We will use this 
example throughout this Section. 

2.2. We complete the set Q with additional variables xj, N < j < N, such that the new set 
Cl forms the smallest algebra, which contains Q and is closed under the Poisson bracket 

{xi,Xj} G fl, for all Xi,Xj G Cl. 

In the above example we add the unit function X2n+i = 1- The complete set contains 
iV = 2?T, + 1 elements satisfying the famous relations 

{XjjXj^n\ {Xj^Yi^Xj^ X2n+1- (1) 
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All other Poisson brackets are zero. 

2.3. We form an iV- dimensional Lie algebra p with the frame {xj}, 1 < j < -/V defined by 
the formal mapping " : x j i— > Xj . The commutators of the frame vectors are formally 
defined by the formula 



We extend the commutator onto the whole algebra by linearity. 

For our example, p is the Lie algebra corresponding to the Heisenberg group (see the 
next Subsection for details). 

2.4. We introduce the algebra ^ of convolutions induced by p. The convolution opera- 
tors are observables in the group quantization, and by analogy with the classical case 
they can be treated as functions of Xj. Particular representations of the convolution 
algebra in spaces lJ{S) give different descriptions of a physical system. The family of 
all one-dimensional representations of ^ corresponds to classical mechanics; various 
noncommutative representations lead to quantum and quantum- classical descriptions 
with different Planck constants. 

For our example the following possibilities exist. 

(a) S = M", Xj — Xj — Mq., Xjj^n — —ihd/dqj, the convolutions are represented 
by pseudo-differential operators (PDO), and we obtain the Dirac-Heisenberg- 

Schrddinger-Weyl quantization by PDO. 

(b) S = R2n^ ^. ^ Xj = Mg., Xj+n — Mp., the convolutions are represented by 
(operators of multiplication by) functions, and we obtain the classical description 
that we started from. 




(2) 



It is an empirical observation that the steps above lead to a nilpotent Lie group, with 
the dual Z of the center Z of the group interpreted as the set of Planck constants. Now 
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we illustrate this fact by a well-known example of quantization, and later in Section ^ by 
constructing a quantum-classical model. 

2.2 The Heisenberg Group Generates Quantum and Classic Me- 
chanics 

In the previous Subsection we claimed that the nth order Heisenberg group H" describes a 
set of quantum particles that constitute an n-degree of freedom system. Here we show how 
this description is achieved. 

M.^ is generated by the n-dimensional translation and multiplication operators e^^'^, e*'''^, 
p,q & ffi" satisfying the Weyl commutator relations 

2-Kip-D 2-Kiq-X 2TTip-q 2mq-X 2mp-D /q\ 

An element of the Heisenberg group g G H" is defined by 2n + 1 real numbers {p,q,s), 
p,q & M", s G M. The composition of two elements g and g' is given by 

{p,q,s){p',q',s') = {p + p' ,q + q' , s + s' + ^{pq' - p'q)). 

Dj, Xj, and / form a 2n + 1 dimensional basis of the Heisenberg algebra [)" with a one- 
dimensional center Z = {si; s G M}. Since all second and higher order commutators of the 
basis elements vanish, H" and f)" are step two nilpotent Lie group and algebra respectively. 

The unitary irreducible representations of the Heisenberg group are classified by the 
Stone- von Neumann theorem [^. They are parameterized by a real number h, the character 



of the one-dimensional center Z. A non-zero h gives non-commutative unitary representa- 
tions acting on the Hilbert space L^(]R") 

PhMP^q,s) = e2-(^-^^+^-^+-^^). (4) 

The n components of X and hD are the usual quantum mechanical position Xj (multipli- 
cation by Xj) and momentum hDj {—ih times differentiation with respect to Xj) operators 
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characterized by the Heisenberg commutator relation 

[hD,,Xk] = -ih5jkl. (5) 

In the limit of zero h the center Z of the Heisenberg group vanishes, and H"' becomes 
isomorphic to M^". The irreducible representations of the latter are homomorphisms from 
M^"' into the circle group acting on C 

Ph=o(p,g) = e2-(P'=+''^). (6) 

The dual H" as a set is equal to {R \ 0} U M^" (see Figure m). It has the natural topology 
coinciding on {M \ 0} with the Euclidean topology. Any sequence of representations {phj}, 
hj 0, hj 7^ is dense in whole M^". The last property is fundamental for the correspondence 
principle. 

The unitary representations of H" can be extended to the convolution algebra L^(H"). 
Namely, if A G L^(EI"), then it defines a convolution on the Heisenberg group 

A -1(9)= [ A{g')b{g^g')dg'. 

The representation ph maps the convolution to the operator 

Ph{A) = I A{g)ph{g)dg 



j j j MP,Q,s)phiP,q,s)dpdqds. (7) 



The p-mechanical equation of motion (see ||20| for details) for an element A(g) (g = {p, q}) 
of the convolution algebra is defined by 

dA{g) 



dt 

with 



2m[H,A]{g) 



[H,A]{g) = [ [H{g')A{g'*g)-A{g')H{g'*g)]dg', 
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The Heisenberg group H 




Figure 1: The Heisenberg group and its dual. 

where H{g) is the Hamiltonian. The non-commutative unitary representations of Eq. (1) 
reduce this equation to the Heisenberg equation of motion for operators acting on the Hilbert 
space L^(]R"). Under the commutative representations of Eq. the p-mechanical equation 
of motion becomes the Hamilton equation for functions on the phase space M^"". 

We consider the last statement in more detail by means of the pseudo-differential calculus 
directly related to these group theoretical developments and the problem of quantization. 
The non-commutative unitary representations of the Heisenberg group allow one to define 
integral operators corresponding to functions on phase space. Given a function a{k, x) on 
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M^" one obtains the operator a{D,X) on L^(R") by the formula 

a{hD,X) = [ J^-'[a]{g)phM9)d9 

= j j J^-^[a]{p,q)e^^'^P-''''+'^-''^dpdq, (9) 

where ^"^[•j denotes the inverse Fourier transform. The trivial integration over s has been 
carried out. The action of the operator o'{D,X) on a function f{x) G L^(R") follows from 
the definition of hD and X [see Eq. (|) and the related paragraph], and is given by 

a{hD,X)f{x) = j J J'^'[a]{p,q)e'''^'"'^^''"'''f{x + hp)dpdq 

= h-^l I J^-^[a]{y^,q)e^^^(^+y^f{y)dydq 

= h-"" I I a{k, ^±y)e^^'^^-y^''/''f{y)dydk (10) 



or 



aihD,X)fix) = / K^ix,y)fiy)dy 



K^{x,y) = j aik,^)e'-'^'^-y^'/''dk, (11) 

where K^r is the kernel of the integral operator a{hD,X). In the language of the pseudo- 
differential calculus the function a{k,x) is called the symbol of the operator a{hD,X). If 
instead of ph^o one uses a commutative representation ph=o, the transformation of Eq. (|) 



reduces to identity and we recover the classical observable a{k,x). Eqs. (||)-(|TT]) are known 
as the Weyl correspondence principle. 

The symbol crtl/iT(k,x) of the product of two operators a^h^ihO, X) = a{hD, X)-T{hD, X) 
can be obtained by application of a non-commutative representation to the convolution on 
the Heisenberg group [see Eq. (^] or directly from the Weyl rule. It is given in terms of the 
symbols of individual operators as 

-2^ 2n 
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It follows from the discussion above that the non- commutative representations of the Heisen- 
berg group transform the p-mechanical equation of motion (|^) into the equation for operators 
on L2(M") 

|a(/^D,X) = ^-^[H,A],ShD,X) (13) 

where [if, Aj^,^ = [i/jj/iA — A'^^H], the operation of taking the product of two symbols jl/i is 
defined by Eq. (p!2D , and the operators A{hD^X) and [if, y4]jj^(/iD, X) are recovered from 
their symbols A{k,x) and [H, A\^^{k,x) by the application of the Weyl transform Eqs. 
(p!l|). This is the quantum- mechanical law of motion in the Heisenberg form. 

In order to obtain the corresponding classical expression it is useful to cast the product 
rule of Eq. (|l2|) in the form of an asymptotic expansion in powers of h. The integration over 
7] and C and the change of variables {u — x) /h u, {v — x)/h ^ v converts Eq. (0) to 



aWr{k,x) = / hF-^[(r]{v,x + uh)J^i[T]{u,x + vh)e^'''^''-''^''dudv, 




where J-'i and J-'^^ denote the Fourier transform and its inverse with respect to the first 
variable only. Expanding a and r in the second variable around x and applying the Fourier 
inversion formula to each term in the Taylor series we obtain 



aUr{k,x) = Yl ' D^kD>{k,x)D'^,D^Mk,x)+0{W 



= J2 —r^^^'-^-'- - Dk,rD^,.YcT{k,x)T{k,x) + 0(/l^), (14) 

where the second subscripts a and r of i^ indicate the symbol to be acted upon. The asymp- 



totic expression for the symbol of the commutator of two operators follows from Eq. (14). 
The even order terms in the sum cancel out to produce 

Hhr - rU(y] {k, x) = 2tJ2 ^"lo' [Dk,aD,,r - i^fc,.i^.,.]'^+ V( A;, x)r (fc, x) 

+ 0{h^^+^). (15) 
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The series expansion of the symbol of the commutator Eq. (15) allows to derive the Poisson 
bracket as the classical limit of the symbol of the Heisenberg commutator of two quantum 
operators 

2ni 

lim—[a^hr -T^h(^]{k,x) = {a{k,x),T{k,x)}. (16) 
/i— >o il 

Since the commutative representations of the Heisenberg group leave symbols of operators 
unchanged, i.e., Jj^n ^'^[(^]{g)ph=o{9)dg = J J J^'^[a]{p,q)e^'''^P^^'^'''>dpdq = a{k,x), we de- 
duce that under the commutative representations the p-mechanical equation of motion (§) 
reduces to the Hamilton equation 

^^A{k,x) = {H{k,x),A{k,x)}. (17) 

In summary, the Heisenberg group contains the exact quantum and classical descriptions 
of a system of particles and provides the correspondence principle between the descriptions. 
We refer the reader to Chapters 1 and 2 of reference for further information on the 
subject. 

3 The Quantum-Classical Equation of Motion 

We proceed to derive an equation of motion for a mixed quantum-classical system. In order 
to do this we look for an abstract mathematical structure that has the same role as the 
Heisenberg group in the case of the standard quantization. The desired structure can be 
constructed based on the following observations. 

First, we need two sets of observables {D,X} and {D',X'} corresponding to quantum 
and classical parts accordingly. We see no reason to assume that an operator from the first 
set does not commute with an operator from the second set. We do assume that each set has 
a Planck constant of its own. Then we let the Planck constant of the second set approach 
zero and obtain the classical limit for the second subsystem leaving the first one quantum. 
We know that "Planck constants" arise as characters of the center. Therefore, we need a 
nilpotent Lie group with a two-dimensional center. 
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This "quantum-classical group" is generated by two sets of variables {hD, X} and {h'D', X'} 
satisfying the commutator relations 

[hDj, Xk] = -ihSjkl, [h'Dj,,Xl,] = -ih'S^kl', l<j,k<n; 1 < /, k' < n' . (18) 

Other commutators are zero. The group has a two-dimensional center Z = {si + s'/'; s, s' G 
M}. Irreducible representations of a nilpotent Lie group are induced by the characters of the 



center p6[. For the quantum-classical group the characters are 

/i : (z, z') exp(i(/iz + h'z')). 

It is clear that for hh' ^ the induced representation coincides with the irreducible repre- 
sentation of EI"+"' on L2(]R"+"'). This corresponds to purely quantum behavior of both sets 
of variables (see Definition |1 . 1| ) . The trivial character h = h' = gives the family of one- 
dimensional representations parameterized by M^Cn+n') ^j^g purely classical description. 
These situations were studied in detail in the previous Section. A new situation appears 
when h ^ and h' = 0, which produce quantum behavior for the first set and classical be- 
havior for the second set. (The h = 0, h' ^ case just permutes the quantum and classical 
parts.) Figure || illustrates these facts. We find that in the topology on the dual to the 
quantum-classical group the quantum descriptions are dense in the quantum-classical and 
classical descriptions, and the quantum-classical descriptions are dense in the classical ones. 



Consider the quantum-classical case in more detail. The quantum-classical representation 
is given by 

PH{p,q,s,p',q',s') = e2-(-'^^+P'^^+''-^+P''='+^'-'). (19) 

where k',x' E M"' and h E 'K \ {0}. In this representation an element of the convolution 
algebra on the quantum-classical group is identified with a quantum-classical operator acting 
of L2(M")®M2"'. The operator can be computed in terms of the Weyl transform of its symbol 



O. V. Prezhdo and V. V. Kisil: Mixing Quantum and Classical Mechanics 14 

Purely quantum Quantum-classical and classical- 

representations quantum representations 




Purely classical representations 

Figure 2: Different types of descriptions generated by a step two nilpotent Lie group with a 
two-dimensional center. 

taken with respect of the quantum (unprimed) coordinates 

a{hD,X,k',x')f{x) = j K^{x,y,k',x')f{y)dy (20) 

K^{x,y,k',x') = h''' j a{k,^^,k',x')e^'''^^-y^''/''dk. (21) 

The quantum-classical analog of the commutator is determined by the limiting procedure 
h ^ 0, h' ^ used to derive the Poisson bracket from the quantum commutator. First we 
need to obtain the expression for the symbol of the product of two operators. We start with 



the expression analogous to Eq. (12), but having two, rather than one set of coordinates. 
Focusing on the primed variables, we carry out the transformations identical to those done 
in deriving Eq. (plGl): 

(TUr{k,x,k',x') = {j^y^ j j j y"rftirf^rfr/dCe^"'[("-")('=-'')-("-")('=-^)]/'^ 

= (^l^"" J J J j dudvdrjdCe^'"'^^''-''^^^-'^^-^'-'-''^^^-^^^/'' 
X (^-^y^ f f dM'dt;'e^-^(^'-''')'=' J-g-i [a] (C, u, v', x' + u'h') 
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X -^3H(??,^^,M ,x +vh) 
•2^ 2n 




X 7^[Dk' ,cr{(;,u,k' ,x')Dx' ,T(ri,v,k' ,x') — D^' ^T{C„u,k' ,x')Dx\c7{ri,v,k' ,x')f 

X [(t(C, u, k', x')T{ri, V, k', x')] + 0{h'^). (22) 



At this point it is safe to drop the explicit dependence of the symbols on the primed variables, 
since the non-commuting nature of the symbols is accounted for by the unprimed variables 
alone. By taking the limit /i' ^ we obtain the expression for the symbol of the quantum- 
classical commutator 



or 



X WiC, u)T{ri, v) - T{C,u)a{r],v) 

+ fe/^^ (2j + l)! i^'='.'^(C-«)^-'.-(^--) 

- /^fc',.(c,.)/^x',.(,,.)]^'^'+'^[a(C,w)r(r7,t;)-r(C,M)a(r/,t;)]} (23) 



[o-,r]s,(fc,x) = (^^y^ J J J J dudvdr^dCe'-'^^^-''^^''-^^-^^^^^^ 

X criC,u)T{ri,v) - T{C,u)a{r],v) 

ih / da{C,u)dT{r],v) drjCu) da{ri,v) \i 

271 \ dk' dx' dk' dx' )\' ^ ^ 

Eq. ( P^ for the symbol of the commutator, together with the rule for calculating operators 
from their symbols given by Eqs. ( ppj ) and (pi]), leads to the following equation of motion 
for a mixed quantum-classical system 

d 27r? 

-A{hD,X,k\x') = —[H,A]i,ShD,X,k',x'). (25) 



O. V. Prezhdo and V. V. Kisil: Mixing Quantum and Classical Mechanics 



16 



This formula determines the evolution of operator A, which depends on quantum and classical 
position and momentum variables and acts on L^(M") (g) R^". The kernel of the operator 
with respect to the L^(R") subspace is given by Eq. (^). 

The quantum-classical equation of motion (|25|) exhibits many desired features. If A and H 
depend solely on quantum or classical variables, Eq. (^) reduces to the corresponding purely 
quantum and purely classical equations, Eqs. (|13D and ([T7|) respectively. Since the quantum- 
classical bracket of the right side of Eq. p5| ) was obtained by selecting a representation for 
the Lie bracket of a Lie group, it has all the properties of Lie brackets, in particular, it is 
antisymmetric and satisfies the Jacobi identity. It is antisymmetric because its symbol is an- 
tisymmetric: the integrand of Eq. changes sign under the permutation a ^ t. To prove 
that the quantum-classical bracket satisfies the Jacobi identity we consider the expression for 
the symbol of the product of three operators a{hD,X, k',x')T{hD,X, k',x')(l){hD,X, k',x') 
By successive application of Eq. ( p2D we get 

^ dudvdrid([a{(, u)T{ri, v) — t{(, u)a{r], v)] 




h' 

_ Q4:TTi[{r}i-u)(vi-r])-{r]i-v){vi-Q]/h^^^^^ ^^-jj 

+ inh'{[Dk'a{C, u)T{r], v) + a{C, u)DyT{ri, v) 

- Dk,T{C, u)a{r], v) - r(C, u)Dk,a{ri, v)] 

X [D^'Cr(C, u)T{r], v) + (t(C, u)D^>T{r], v) 

- D,,r(C, u)a{7], v) - r(C, u)DMv, v)]} + 0{h'^)}. (26) 

After taking the limit lim/i/^o jp- of the first and higher order terms in h' and performing the 
integrations, this cumbersome expression can be rewritten in a more compact symbolic form 
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as 



crjjrtl^ = GT(j) - (j)To + i{ak'T(j)x' + OTk'cjyx' - rkiocj)^' - TOk'(t)x' 

-(pk'CTx'T - 4>k'Crrx' + (pk'Tx'Cr + 4>k'T(Tx'), (27) 

where the subscripts k' and x' indicate differentiation with respect to these variables, and 



the ordering of the symbols is to be kept track of. Given Eq. (^7]) it is straightforward to 
check that the Jacobi identity holds for the symbols of operators: [[a, r]|j^, 0]jj,^(fc, x, A;', x'), 
and, therefore, for the quantum-classical bracket. 

We used the Weyl correspondence principle to link quantum and classical mechanics and 
to derive the quantum-classical equation of motion. It is well known, however, that the Weyl 
correspondence in not unique in mapping phase space functions to Hilbert space operators. 
In fact, there exist arbitrary many such mappings, differing in the order assigned to products 
of position and momentum operators. The Weyl rule leads to symmetrized products. For 
example, it maps kx to ^{hD ■ X + X ■ hD) . Another variant of the correspondence principle, 
widely used in the mathematics community because of its simpler form, is due to Kohn 
and Nirenberg. It keeps momentum operators on the right, mapping kx to X ■ hD. It is 
straightforward to obtain a mixed quantum-classical equation of motion using the Kohn- 
Nirenberg calculus. One starts with ph = ^'^'^n-x ^2mp hD instead of ph = e^'^^^'^^'^'^'^^ and 
follows the same steps. The result is 

^^A{hD,X,k',x') = -^[H,A]l''ihD,X,k',x') (28) 

with the correspondence rule 

a{hD,X,k',xT''fix) = J K^''ix,y,k\x')fiy)dy, 

A'f^(x,i/,A;',x') = /"a(A;,x,fc',x')e2"*(^-^)'=/'^dA; (29) 



and the following formula for the quantum-classical bracket of two symbols 

■2^ 2n 



]^^^{k,x,k',x') = (-) dudve^^'^''-''^^''-^^'^'^a{v,x)T{k,u)-T{v,x)a{k,u) 
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ih 



da{v,x) dT{k,u) driv^x) da{k,u) 



27rL dk' dx' dk' dx' J J' ^^^^ 

These expressions are somewhat simpler than those obtained by the Weyl correspondence. 
Eqs. (^,^) and (P^P^D are preferable, however, as they preserve the simplectic invariance of 
the phase space variables and lead via the Wigner transform from the density matrix to the 
quantum quasi-probability function that is closest to the classical probability density [p^ . 



4 Quantum-Classical Coupling for Harmonic Oscilla- 
tors 

The standard quantization can also be obtained by application of a projection to the classical 
system under consideration, which leads to a simple description of harmonic oscillators. We 
briefly summarize this topic and apply it to the case of coupled quantum and classical 



oscillators. Further information can be found in |28, 2£, 32, 33, 3|| and references 

therein. 

Let L2(C",(i/in) be a space of functions on C" square-integrable with respect to the 
Gaussian measure 

where dv{z) = dxdy is the Euclidean volume measure on C" = M^". The Segal-Bargmann |^ 



^ (or Fock) space F2(C") is the subspace of L2(C", c?/i„) consisting of all entire functions, 
i.e., functions f{z) that satisfy 

1^ = 0, 1 < J < n. 



Denote by Pq the orthogonal Bargmann projection |^ of L2(C"', onto the Fock space 
F2(C"). Then 

k{q,p) Tk(q+ip) = Pgkiq + ip)I (31) 

defines the Berezin-Toeplitz (anti-Wick) quantization, which maps a function k{q,p) = k{q + 
ip) on M^*^ = C" to the Toeplitz operator T^. There exists an identification between Berezin 
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and Weyl quantizations 31, 3^]. The identification has an especially transparent form 



for the observables depending only on p or q. 

The Berezin-Toeplitz quantization is related to the Heisenberg group more intuitively 
than the representation of Eq. (^). On a geometrical level, consider the group of Euclidean 
shifts a : 2; I— > 2; + a of C". To obtain unitary operators on L2(C", dfi) the shifts are multiplied 
by the weight function 

a : f{z) ^ f{z + a)e-^™/2_ ^32) 

This mapping determines |^ a unitary representation of the (2n+l)-dimensional Heisenberg 
group acting on ^2(0", dfi). The mapping preserves the Fock space F2(C"'), and hence all op- 
erators of the Eq. (^) form commute with Pq. The operators are generated by infinitesimal 
displacements 

n 



Ell/ I ■ ii\ I II I ^ II , ■ n 

S - - '^3 ) + «i - ^3 + 

k=i ^ ^ i ^ i 



where a = (ai, . . . , a„), z = {zi,...,Zn) G C" and aj = {a'j,a"),z = {z'-,z'-) G M^. The 
generators form a linear space with the basis 

4' =7 fir--;--?). 4"=7fj|r--;+'-i')- (33) 



i \ dz'j ■' ^ J ^ i \ dz" 



The basis vectors commute with the Bargmann projector Pq. The operators 



4' = 7(|r-4 + -;). 4" = 7(^-<-m (34) 



commute with the basis vectors, and we anticipate that Pq produces a self-adjoint represen- 
tation of convolutions with respect to Xj', ^j", and unit operators. Indeed, 

Proposition 3 ||^ The Bargmann projector Pq defines a representation of convolutions 
induced by the Weyl-Heisenberg Lie algebra f)„ operating on by Eqs. (0). The kernel 
b{t, C), t G M, zeta G C" of the representation is given by the formula 

6(t,C)=2"+i/2e-ie-(*'+^^"/2). 
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We move on to apply the Bargmann projection technique to the quantum-classical cou- 
pling of harmonic oscillators. 

Example 4 In the Segal-Bargmann representation the operators of creation and an- 



nihilation of the jth state of a particle are = Zjl and a^- = d/dzj, correspondingly. 
Consider an n degree of freedom harmonic oscillator with the classical Hamilton function 



1 " 



2 

The corresponding quantum Hamiltonian is obtained by the Bargmann projection 

\ ^ \ ^ d 

THi,,p) = ^Pq + P?)/ = 2 K + E ^^■^)- (35) 

j=i j=i J 

The right side of Eq. (^) is the celebrated Euler operator. It generates the well known 
dynamical group Chap. 1, Eq. (6.35)] 

e^tT,,,,,) j(^) = e^"*/2/(e^*z), f{z) G F2, (36) 

which induces rotation of the C" space. The evolution of the classical oscillator is also given 
by a rotation, that of the phase space M^"' 

z(t) = GtZo = e**2;o, z{t) = p{t) + iq{t), zq=pq + iqo. (37) 

The projection Pq leads to the Segal-Bargmann representation, providing a very straight- 
forward correspondence between quantum and classical mechanics of oscillators, in contrast 
to the rather complicated case of the Heisenberg representation ||38|, Prop. 7.1 of Chap. 1]. 

The powers of z are the eigenfunctions 0„(-2) = z"- of the Hamiltonian (^), and the 
integers n are the eigenvalues. Either pure or mixed, any initial state of the oscillator 
remains unchanged during the Eq. ( ^61) evolutions and no transitions are observed. 

Now consider classical and quantum oscillators coupled by a quadratic term 

H{p,q;p',q') = ^{p'^ + p^ + x'^ + + a{x' - xf). (38) 
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Applying the canonical transformation (see [3£, § 23. D]) 

/ ?i + 92 qi-q2 , Pi+ P2 Pi - P2 .on^ 

Q = — , 9 = — , P= — p — , P= — F — , (39) 
V2 V2 V2 V2 

or, equivalently, introducing complex variables z = q + ip, z' = q' + ip', Zi = qi + ipi, 

Z2 = q2+ ip2 

/ _ Zi + Z2 ^Z\ — Z2 

z' + z z' — z 

we get rid of the coupling term 

H{pi,qi]P2,q2) = ^{pl+pl + i^lql + ^lql), (41) 



where tui = 1, =vT+2a. The two uncoupled oscillators evolve independently 

z,{t) = e'^-^'zM, Z2{t) = e'^-^'zM 

The dynamics in the original coordinates, however, is not trivial. The primed and unprimed 
(quantum and classical) variables mix 

^ \'^) ~ 2 ' ( 

z[t) — - . [AS) 

Suppose that the classical subsystem is initially localized at a point Zq the phase space 
and the quantum subsystem is in its n-th pure state 0(z', z; to) = 6{zq — z') (g) z"-. The 
dynamics of the combined system is given by 

(P{z',z;t) = 6izo-^ ^-j (44) 

o2iuiit g2ia)2t^^' _|_ ^g2ja;it _j_ ^2iijJ2t'j^^ 
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During the evolution the classical subsystem, [Eq. (^^1)] is always sharply supported, i.e. 
represented by the delta function, while the quantum subsystem, [Eq. (^5D] evolves into 
a mixed state. The binomial ( ^5]) contains all powers of z less or equal to n (^z^ < n). 
Therefore, there exists a non-zero probability for the quantum subsystem to make a transition 
from the initial pure state into a lower energy state , k < n. It is remarkable that in 
this particular case the interaction with the classical subsystem can only decrease the initial 
energy of the quantum one. The overall dynamics is (quasi)-periodic with the recurrence 
time determined by the frequencies uji and uj2. 

5 Discussion and Conclusions 



The quantum-classical equation of motion (|25| ) can be applied in several ways depending on 
how one describes the quantum and classical subsystems. The classical subsystem can be 
treated on the level of trajectories, in which case it is represented by a point in the phase space 
{k[,x'j} evolving according to the Hamilton equations with the quantum-classical bracket of 
Eq. (^) regarded as a modification of the Poisson bracket. If at the same time the Heisenberg 
equation of motion with the quantum-classical bracket in place of the commutator is used 
to describe the evolution of quantum operators, the mean field approximation is recovered. 
Namely, the quantum mechanical average of Eq. (pS]) with respect to the wave function \l/ 
is given by 

(9 27r? 

^^mA{k\x'm = —{nH,AUk',x')\^). m 

If A is a purely quantum mechanical observable independent of classical variables, the deriva- 
tives dA/dk' and dA/dx' in the quantum-classical bracket Eq. (^) vanish, and we obtain 

-(vi/|A|vl/) = -^{nH{k',x'),A]\^) (47) 

with the Hamiltonian H parametrically depending on the classical phase space variables k' 
and x' . Substituting the variables in place of A{k'^x') in Eq. (|46|) , we recover the classical 
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equations of motion, the classical Hamiltonian being the quantum mechanical average of the 
total Hamiltonian 

dk' d{-^\H{k\x')\^) 



dt dx' 
dx' d{^\H{k',x')\^) 



(4J 



dt dk' 

Eqs. (|47|) and (^81) constitute the traditional mean field approximation: classical variables 
are coupled to the expectation values of quantum observables. 

The quantum-classical equation of motion can be looked upon in a different way, namely, 
as a Liouville-von Neumann equation for a mixed distribution p{hD,X,k',x'). Selecting 
a quantum basis we get a set of coupled equations for classical phase space distribution 
functions pij{k',x') corresponding to each pair of the quantum basis states i, j 
dpij _ 27rz^r ihrdHl^dpkj dp*^dHkj\l 



k 



In the purely quantum and classical limits the two derivatives [Eqs. (^), ( ^8|) and 
Eq. (|49|) ] of the quantum classical equation of motion [Eq. (^)] are equivalent. In the mixed 
case they are not, because of non-local correlations in the classical subsystem induced by 
its interaction with the quantum one. Such correlations, inherent in the Liouville-von Neu- 
mann equation, do not appear on the level of individual trajectories. In particular, if the 
equations (^) are integrated with the initial conditions pu{k', x') = 5{k' — kQ)6{x' — Xq) and 
Pij = 0, V {ij} 7^ {11}, at later times, in general, pij ^ because of the couplings ifjj 7^ 0. 
In other words, classical phase space distribution functions associated with different quan- 
tum states differ and mix. Spreading of initially localized classical distributions is enhanced 
by their mutual mixing. An initially localized phase space distribution pn corresponding 
to the ground quantum state populates excited state distributions. They undergo diverging 
evolutions, and later, when pn gains an admixture of the excited state distributions, it nec- 
essarily becomes delocalized. This, of course, can not happen in the mean field approach, 
where the classical subsystem is described by a single trajectory. Non-local correlations are 
averaged out in the mean field approximation. 
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If there is no coupling between quantum states, phase space distributions neither mix 
nor spread beyond the classical divergence. Consider a purely adiabatic case, where the 
quantum basis states are the instantaneous eigenstates of the quantum Hamiltonian. In 
the absence of non-adiabatic coupling Eq. (|49|) splits into a set of uncoupled equations 
for classical distribution functions corresponding to individual adiabatic quantum states. 
An adiabatic evolution of every distribution function can be equivalently described by the 
classical trajectory mean field approach [Eqs. (P7|)-(^8D] with the wave function \1/ being the 
corresponding adiabatic eigenstate of the quantum Hamiltonian. 

The Liouville-von Neumann equation for the evolution of the mixed distribution reduces 
to the coupled equations of the mean field approximation if the quantum-classical function 
is constrained to 

p = \^){^\.6{k' -k")6{x' -x"), (50) 

where k',k" and x',x" are the n-dimensional classical momentum and position vectors. Under 
this constraint the quantum part of the mixed function always remains a pure state, and 
the classical part is always represented by a delta function. Substituting expression (|50D for 
p in place of A in Eq. (p5D and integrating out the phase space variables (J dk" J dx") we 
obtain the von Neumann equation with the quantum Hamiltonian depending on classical 
coordinates H{k',x'). Since the density matrix entering the equation is constructed from 
a pure quantum state the von Neumann equation for the density matrix is isomorphic 
to the Schrodinger equation and to Eq. (^7|). To recover the mean field equations for the 
classical variables (^8[) we substitute ( ^0]) into (|25|), multiply both sides by either k' or x', 
and integrate over quantum and classical coordinates. 

The mean field approximation can not reproduce the non-local correlations within the 
classical subsystem due to interaction with individual parts of the quantum subsystem. At 
the same time, such correlations naturally appear in the solutions to the Liouville-von Neu- 
mann equation. Unfortunately, the Liouville-von Neumann equation does not provide signif- 
icant computational advantage over the pure quantum von-Neumann equation, since both 
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deal with delocalized distributions. On the other hand, propagation of an individual classical 
trajectory via the Hamilton equations of motion is far less demanding than integration of 
the Liouville equation. The idea of running and then averaging over a few representative tra- 
jectories instead of propagating the total phase space distribution function is fully exploited 
in classical molecular dynamics simulations (see, for instance, reference [^). In order to 
account for the quantally induced non-local correlations among the phase space variables, 
while retaining a trajectory description for the classical subsystem, we interpolate between 
the mean field and Liouville-von Neumann approaches by developing a multiconfiguration 
version of the mean field method. We start with the quantum-classical distribution function 

P = (51) 

i j 

where the sums run over the number of independent configurations. Taking the wave func- 
tions to be orthonormal (\E'j|\E'j) = 6ij we substitute (^) into Eq. (^) and integrate over 
the classical variables (J dk' J dx') to obtain the von Neumann equation for the quantum 
density matrix g = Y.i 

= ~j^'^^[Hikikii.^ Xif^)gik — QkjHkjiki^j, Xf^j)]. (52) 

k 



To get the equations of motion for the classical variables kij, Xij we substitute (^) into 
Eq. (|25|), multiply both sides by the corresponding variable, and integrate over all degrees 
of freedom. We obtain equations of motion for the "diagonal" positions and momenta as 

dt - ^ dx',^ ' 
The corresponding expressions for the "non-diagonal" variables are more complicated 



-Qf = 2^Hi,{hk,Xik)Xkj-XikHkj{kik,Xik) + , (54) 



k 



O. V. Prezhdo and V. V. Kisil: Mixing Quantum and Classical Mechanics 



26 



and similarly for momenta. In order to keep the Hamiltonian matrix || Hij{k^^x^) || her- 
mitian we require that the position and momentum "matrices" are symmetric: k[- = kj^, 
x[j = x'jj^. With this constraint the "non-diagonal" evolutions simplify, the first two terms 
in formula (|5^ disappear, and the dynamics of the "non-diagonal" phase space variables 
coincide with the average dynamics of the "diagonal" variables 

dt ~ 2 \ dt dt 

— ^ \^ f fj* ' — ' n j-^I^itj-U* ' — ' n I ^-^jfc 



2 V ^ ^''ik ^^k ^ 2\ dt dt 

~dt ~ 2\dt^~dr)- ^^^^ 

Apart from assigning the "non-diagonal" coordinates and momenta unique values a simple 
physical meaning: x'^ = (x-j + x'jj)/2, k'^j = {k[^ + kjj)/2, Eqs. ( ^5|) reduce the number of 
independent classical trajectories in the n-dimensional multiconfiguration mean field approx- 
imation from to n. 

The idea of introducing multiple configurations is not new. Apparently, it originated 
in quantum chemistry as an improvement on the self-consistent field solution of the time 
independent Schrodinger equation. The time dependent fully quantum multiconfiguration 
self-consistent field approach is discussed in reference Equations of motion similar 

to our version of the quantum-classical multiconfiguration approximation were proposed by 
Diestler 0]. The method was devised to account for non-local correlations among the classical 
degrees of freedom of a condensed phase solvent interacting with a quantum solute. Our 
approach has a different derivation and is less computationally demanding by a factor of n. 
The number of configurations in the multiconfiguration mean-field method does not have 
to be the same as the number of quantum (adiabatic) basis states. The former is usually 
less than the latter, and, in the case when the classical phase space can be separated into a 
union of several weakly connected regions, it should be determined by the number of such 
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regions. For example, a double well system would require two configurations — two classical 
trajectories, each originating in its own well. 

The quantum-classical equation of motion is derived here via the simplest extension of 
the Heisenberg group of quantum and classical mechanics to the group allowing for two sets 
of variables, such that the variables of the first set do not commute, while the variables of 
the second set do. Obviously, there exist other groups satisfying this requirement. It is 
possible, for instance, to consider a step three nilpotent Lie group and the corresponding 
algebra decomposable as a vector space into the three subspaces Vq, Vi, and V2 having the 
following properties: The elements of Vq commute with all elements and form the center. 
Commutators of vectors from Vi belong to Vq. Commutators of vectors from V2 belong to Vi. 
By taking a representation of this group that maps the center Vq to zero, we would obtain 
another model for a mixed quantum-classical system, where vectors from Vi would corre- 
spond to classical degrees of freedom, since their commutators vanish, while vectors from V2 
would describe quantum variables. It is likely, though, that this scheme will exhibit proper- 
ties atypical for quantum and classical mechanics, since step three nilpotent Lie groups differ 
from step two groups and the Heisenberg group in particular. There exist, however, some 
advantages in dealing with general nilpotent Lie groups. For example, the relativistic quan- 
tization considered in reference |^ is based on a representation of the simplest step three 
nilpotent Lie group [meta Heisenberg group |^) spanned by the Schrodinger representation 
of the Heisenberg group and the operators of multiplication by functions. Application of the 
quantization rules to the appropriate Lie algebras leads to quantum-classical constructions 



for string theory, conformal field theory, and Yang-Mills theories ^4 



In summary, we presented a systematic approach to coupling quantum and classical de- 
grees of freedom based on a generalization of the unified description of quantum and classical 
mechanics in terms of convolutions on the Heisenberg group. Considering the simplest ex- 
tension of the Heisenberg group that allows for two distinct sets of variables, we derived a 
quantum-classical equation of motion. The quantum-classical bracket entering the equation 
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is a Lie bracket and, therefore, possesses the two most important properties common to the 
quantum commutator and the Poisson bracket: It is antisymmetric and satisfies the Jacobi 
identity. We explicitly constructed the quantum-classical equation of motion for coupled har- 
monic oscillators and discussed approximations to the equations applicable to more general 
cases. 
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